Structural disorder, magnetism, and electrical and 
thermoelectric properties of pyrochlore Nd2Ru207 



Michael W Gaultois^ Phillip T Barton^ Christina S Birkel^ 
Lauren M Misch^, Efrain E Rodriguez^, Galen D Stucky^ 
and Ram Seshadri^ 

^ Materials Department, Department of Chemistry and Biochemistry 

and Materials Research Laboratory University of California, Santa Barbara, CA 

93106, USA 

^ Department of Chemistry and Biochemistry 
University of Maryland, College Park, MD 20742, USA 

E-mail: mgaultois@mrl.ucsb.edu 

Abstract. Polycrystalline Nd2Ru207 samples have been prepared and examined using 
a combination of structural, magnetic, and electrical and thermal transport studies. 
Analysis of synchrotron X-ray and neutron diffraction patterns suggests some site 
disorder on the A-site in the pjnrochlore sublattice: Ru substitutes on the Nd-site 
up to 7.0(3)%, regardless of the different preparative conditions explored. Intrinsic 
magnetic and electrical transport properties have been measured. Ru 4d spins order 
antiferromagnetically at 143 K as seen both in susceptibility and specific heat, and 
there is a corresponding change in the electrical resistivity behaviour A second 
antiferromagnetic ordering transition seen below lOK is attributed to ordering of Nd 
4f spins. Nd2Ru207 is an electrical insulator, and this behaviour is believed to be 
independent of the Ru-antisite disorder on the Nd site. The electrical properties of 
Nd2Ru207 are presented in the light of data published on all A2RU2O7 pyrochlores, 
and we emphasize the special structural role that Bi^+ ions on the A-site play in 
driving metallic behaviour High-temperature thermoelectric properties have also 
been measured. When considered in the context of known thermoelectric materials 
with useful figures-of-merit, it is clear that Nd2Ru207 has excessively high electrical 
resistivity which prevents it from being an effective thermoelectric. A method for 
screening candidate thermoelectrics is suggested. 



PACS numbers: 75.47.Lx 75.50.Ee 75.50.Lk 72.20.Pa 
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Figure 1. Interpenetrating a (a) "cubic-ice" diamondoid A2O' lattice with (b) a lattice 
of corner-connected BOq octahedra gives the (c) ^12-82060' pyrochlore structure of 
Nd2Ru207. Space group FdSm [No. 227, A at (i, i, i), B at (0,0,0), O at (xq, g, |), 
and O' at (|, |, |)]. Black spheres are Nd, blue octahedra surround Ru, and the orange 
spheres are O. 

1. Introduction 

The A2-B2O6O' pyrochlore structure, shown in figure[T| comprises two interpenetrating 
B2OQ and A2O' sublattices lH]]. The electrical properties of pyrochlore ruthenium 
oxides (y42Ru207, A= Pr through Lu, Y, and Bi) are of long-standing interest. In 
y42Ru207 ruthenates, electrical conductivity is expected to take place through the 
RU2O6 network of RuOe octahedra, ie. the -B2O6 sublattice. While all the rare-earth 
members are insulating, Bi2Ru207 is metallic. Several investigators have sought to 
explain the difference in conductivity by examining changes in Ru-O-Ru bond angle 
— determined by the combination of Ru-0 bond length and the size of A — and 
consequent changes in orbital overlap and bandwidth ^ |3|] . Electron spectroscopic 
investigations of Y2RU2O7 compared with Bi2Ru207 led Cox et al. flm to conclude 
that the principle difference seen is the participation in conduction of Bi 6s states 
in Bi2Ru207. This participation broadens the Ru 4d band width sufficiently that a 
metallic ground state is preferred over one that is correlated and insulating. However, 
Shoemaker et al.[5] contrasted the computed electronic structures of insulating 
Bi2Ti207 with that of conducting Bi2Ru207 and found no difference in the presence 
or absence of Bi 6s states near the Fermi energy in these two compounds. 

Pyrochlore-type rare-earth ruthenium oxides are also of great interest for their 
magnetic properties. The sub-lattices of corner-connected A4 and B4 tetrahedra can 
result in magnetic frustration when either A or B is separately magnetic ^ . While 
magnetic properties of Nd2Ru207 have been reported previously |]3, we show here 
that the presence of ferromagnetic impurity phases compromised the earlier analysis. 
There has also been some confusion in the literature about the nature of the ordering 
of the Ru 4d spins in Nd2Ru207; several reports have suggested that the ordering 
is glassy [H [H [lOl [IH [l2ll • We find compelling evidence that Nd2Ru207 is not a 
spin-glass, and propose a reason for the observed history-dependence of the ZFC-FC 
magnetic susceptibility and weak ferromagnetism below the ordering temperature. 
Heat capacity studies of the magnetic transition in Nd2Ru207 are also presented here. 

In A2RU2O7 pyrochlores the structural modifications between metallic and 
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semiconducting members are small, suggesting that the electrical transport properties 
are delicately positioned near the edge of localized and itinerant behaviour. This 
positioning, at the edge of metal-insulator divide, potentially serves as a prime 
locator for thermoelectric materials: a high Seebeck coefficient is generally found 
in insulators, with potentially acceptable electrical conductivity on the metallic side, 
striking the right balance in properties [lT3l [T4ll . For this reason, we investigate 
the high-temperature thermoelectric properties of Nd2Ru207. When compared to 
well-known chalcogenide thermoelectric materials, for e.g., PbTe, Bi2Se3, and Bi2Te3, 
metal oxides are of interest due to the expectation of higher stability at elevated 
temperatures and the prospect of using less toxic and more abundant elements. The 
discovery of high thermoelectric performance in Na2;Co04[15] has led to renewed 
interest in oxide thermoelectric materials, despite the figure of merit (ZT) of this 
system and other bulk oxide materials being too small for widespread use ||T4l [T6ll . 
Towards the goal of a more directed and effective search for high-performance 
metal-oxide thermoelectric materials, we also introduce a new type of plot for data 
visualization that can be used to rapidly screen potential candidates for thermoelectric 
performance. 

2. Experimental details 

Polycrystalline samples were made by direct reaction of constituent oxide powders 
(RUO2, 99.99%, Sigma-Aldrich; Nd203, 99.99%, Alfa Aesar). Owing to the volatility 
of Ru oxides in air at high temperatures (>1040°C), samples were prepared with a 
1 mol% excess RUO2. Pellets were cold-pressed and annealed at 1000°C for 7 days 
with several intermediate grindings, after which starting materials were still found by 
lab XRD. Following this initial reaction, pellets were wrapped in Pt-foil and annealed 
in evacuated silica ampoules at 1100°C for 7 days. The samples were cooled to 800°C 
at 0.5°C/min, and further annealed at 800°C for 7 days to promote healing of defects. 
Finally, samples were cooled slowly to room temperature at 0.5°C/min. If Nd2Ru207 
is annealed in air at 1040°C or 1060°C for 12 days, Rietveld analysis of XRD patterns 
indicates the sample contains 78 mol% NdsRuOy and 22 mol% Nd2Ru207 when the 
annealing temperature is 1040°C, and pure NdaRuOy when the annealing temperature 
is 1060°C. 

To avoid the formation of NdsRuOy at the high temperatures and extended 
annealing times required for solid state reactions, polycrystalline samples of >99% 
purity were also made by ultrasonic spray pyrolysis (USP) of aqueous nitrate 
precursors [Ru(NO)(N03)3, Sigma-Aldrich; and an aqueous solution of Nd^+, 
prepared by dissolving Nd203 in ^2M HNQcil UTTl \T8L A modified ultrasonic 
humidifier was used to generate a mist, which was passed through a tube furnace 
(30.5 cm heating zone) at 700°C with a positive pressure of air at a flow-rate of 5 scfm 
(;^140L/min). The product was collected in water, and the mixture was evaporated 
at 70°C; the resulting black solid was cold-pressed into a pellet and annealed in air at 
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1100°C for 8 hours. Samples were then quenched in air (removed from the furnace 
at 1100°C) or annealed in air for 7 days at 800°C then cooled to room temperature 
at l°C/min. Rietveld refinement of the structure with synchrotron XRD data revealed 
no significant differences in the bulk long-range structure between samples that were 
quenched rapidly and samples that were annealed at intermediate temperatures and 
cooled slowly, nor were there differences in structure between samples made by USP 
and samples made by solid state reaction. USP offers a rapid preparatory route to 
prepare Nd2Ru207 in ^24 hours, rather than ^14 days. 

Laboratory XRD was performed using a Philips X'Pert diffractometer with Cu Ka 
radiation and using mis-cut Si sample holders to reduce background signal. High- 
resolution synchrotron X-ray diffraction (XRD) data on finely-ground powder was 
acquired at 100 K and 295 K at beamline 11-BM at the Advanced Photon Source 
(APS), Argonne National Laboratory, using an average wavelength of 0.413 A on a 
diffractometer that has been described in detail by Wang et al. [fT9i] . The precise 
wavelength was determined using a mixture of Si (SRM 640c) and AI2O3 (SRM 
676) NIST standards run in a separate measurement. Samples were contained in 
0.4 mm diameter Kapton capillaries and the packing density was low enough that 
absorption was not noticeable. Neutron powder diffraction was performed on the BT- 
1 diffractometer at the NIST Center for Neutron Research. A Cu(311) monochromator 
was used, with a constant wavelength of A = 1.5402(2) A and a second-order 
contribution at A/2. Data was collected at 300 K over the range of 3° to 168° 26 
with a step size of 0.05°. All diffraction data shown here are from a sample made by 
solid state reaction, though Rietveld analysis of data obtained by synchrotron XRD on 
other samples made by solid state reaction or USP lead to identical results. 

X-ray total scattering was performed at beamline 11-ID-B at the APS, using a 
wavelength of 0.137020 A. The pair distribution function (PDF) was extracted with 
PDFgetX2 [l20ll using Qmax = 28 A \ and full-profile PDF structure refinement was 
completed using PDFgui (2T]. i?wf and xled were refined on 1950 data points. 
Key instrumental parameters were Qbroad = 0.0551 A"^ and Qdamp = 0.00963 A"^ 
respectively, using a Ce02 standard run in a separate measurement. 

Magnetic properties of powders were measured using a Quantum Design MPMS 
XL-5 SQUID magnetometer. In addition to DC measurements, frequency-dependent 
AC measurements were performed in a small temperature range between 130 K 
and 155 K. Low-temperature electrical transport properties and heat capacity were 
measured using a Quantum Design Physical Properties Measurement System. Samples 
for electrical transport measurements employed the 4-probe geometry on a pellet of 
sintered powder with dimensions of approximately 9mmx 3 mm x 3 mm. Electrical 
contacts were made with copper wire and silver epoxy. Three samples were run 
to ensure reproducibility: two made using USP with apparent densities of 38% 
and 53%, and one made by ceramic preparation with an apparent density of 58%. 
Electrical resistivity at 300 K varied from 1.2ficm to 2.4^2 cm between the samples, 
and is consistent with previous reports of 1.8f2cm [l22ll . Despite slight differences 
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in magnitude of resistivity, all samples displayed consistent temperature-dependent 
behaviour. Heat capacity measurements were collected on a pellet of mass 13.90 mg 
made by pressing 50wt% of sample and 50wt% of silver powder (99.99%, Sigma- 
Aldrich) and analyzed using the thermal relaxation dual-slope method. A thin layer 
of Apiezon N grease ensured thermal contact between the platform and the sample. 
The heat capacity of the Apiezon N grease and silver were collected separately and 
subtracted from the measured sample heat capacity. 

High-temperature thermoelectric properties (electrical resistivity and Seebeck 
coefficient) were measured with an ULVAC Technologies ZEM-3. Sample pellets 
had approximate dimensions of 9 mm x 3 mm x 3 mm. Measurements were performed 
under a helium under-pressure, and data was collected through three heating and 
cooling cycles to ensure sample stability and reproducibility. No changes in physical 
properties were observed between cycles, and analysis of the lab XRD pattern of the 
materials after measurements showed no changes in structure nor new phases. Two 
different samples (USP and ceramic) were tested to verify consistency. 

3. Results and Discussion 

3.1. Structure 

The high symmetry ideal pyrochlore crystal structure is completely determined by 
the cubic cell parameter and a single positional structural parameter, associated by 
the position x of O which is sited at (xo, |, |)- However, there are many types of 
disorder that arise in this structure type, and a careful examination of the structure is 
necessary to understand properties. For example, Vanderah et al. recently described 
the widespread presence of antisite substitution on the A-site; up to 25% of the large 
A-sites can be replaced with small B-site metal ions [l23ll . Moreover, local off-centering 
of A-site cations is well-known in Ru pyrochlores such as Bi2Ru207 p5l, 124']. The 
presence of disordered oxygen vacancies has been observed in metallic members of 
A2_2:Bix.Ru207 solid solutions [25]. Due to the relatively small X-ray scattering factor 
of O compared to the other elements present in y42Ru207, use of only lab XRD has led 
to inaccurate determination of the O atomic position, as demonstrated by Kennedy 
and Vogt [f25ll . This inaccuracy has a large effect on the reported Ru-O-Ru bond 
angles, which are known to critically influence electrical properties. For example, in 
Bi2Ru207, there is a 6° discrepancy, and for other A2RU2O7 members there is up to a 
2° discrepancy in the reported Ru-O-Ru bond angles []3l |24l [26ll . 

Combined Rietveld refinement was carried out using room-temperature data 
sets (figure|2]), where the structure was refined on 51295 data points. During 
the combined refinement, the synchrotron X-ray wavelength was fixed while the 
neutron wavelength was allowed to vary, though the refined wavelength was within 
two standard deviations of the starting value determined by previous instrumental 
calibration. Isotropic displacement parameters were used to describe the electron 
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Figure 2. Combined Rietveld refinement of the structure with (a) room-temperature 
synchrotron XRD with A = 0.41295A and (b) neutron diffraction with A = 1.5402A. 
From the refinement, the cell parameter a is determined to be 10.342312(8) A at room 
temperature, and the position xq corresponding to the O atom is 0.33012(7). 



density of the atoms in the structure, as refinement of anisotropic displacement 
parameters did not improve the quaUty of the fit. Rietveld refinement was performed 
with XND code [|27ll , and structures were visualized using VESTA [l28ll . 

Rietveld refinement of the structure using synchrotron powder XRD (figure|2](a)) 
data indicates samples were 99.5(1) mol% Nd2Ru207 and 0.5(1) mol% RUO2. The 
minor presence of RUO2 is not expected to influence the physical properties reported 
here. There were no significant differences in structure between samples that were 
quenched rapidly and samples that were annealed at intermediate temperatures 
and cooled slowly. Additionally, there were no significant differences in structure 
between samples made by USP and samples made by solid state reaction, nor 
were there differences between samples annealed under static vacuum and samples 
annealed in air. Taken together, the long-range structural order of Nd2Ru207 appears 
to be insensitive to preparation conditions and methods. In particular, annealing 
Nd2Ru207 under low oxygen partial pressures (i.e., in evacuated ampoules) does 
not appear to lead to oxygen deficiency. More evidence of this stability is presented 
in a later subsection, as electrical resistivity does not change after several heating 
cycles between 300 K and 900 K under oxygen-free conditions. This suggests that O 
vacancies not already present in the ordered pyrochlore structure are not formed to an 
appreciable extent in Nd2Ru207. To check for the presence of oxygen vacancies, the 
site occupancy of O' was allowed to refine freely; the best refinements were consistent 
with complete occupancy. This finding is supported by previous neutron-diffraction 
studies, which found no evidence for additional oxygen vacancies in Nd2Ru207 [i25i1 . 
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Figure 3. Refinement with A-site disorder converges to 7.0(3)% Ru on the A-site and 
sHghtly improves the fit to the data, as demonstrated above by the better description of 
the peak-shape. Changes in the calculated diffraction pattern due to site disorder are 
small, with the most significant and diagnostic change occurring at the 311 reflection. 
For comparison, the most intense peak in this pattern has 85 000 counts. 



The oxygen positional parameter {xq, |, |) converged to xq = 0.33012(7), identical 
to previous reports using neutron diffraction [l25l [29ll . 

Synchrotron X-ray diffraction studies offer the advantage that the Nd^+ and 
Ru^+ X-ray scattering factors are sufficiently distinct due to the large difference 
in atomic numbers. In contrast, neutron diffraction does not clearly distinguish 
between the similar coherent scattering lengths of Nd (7.03 fm) and Ru (7.69 fm) 
[l30ll . The synchrotron diffraction data allowed the recent findings of Vanderah et 
al. regarding antisite disorder on the cation sites to be tested. The site occupancy 
and atomic displacement parameters are often correlated, so synchrotron X-ray 
data was collected at both 100 K and room temperature (295 K). However, these 
parameters were not strongly correlated in this investigation, and a combined Rietveld 
refinement using multiple temperatures did not change the outcome of the analysis. 
Accordingly, we continue our discussion using the combined Rietveld refinement of 
room-temperature X-ray and neutron diffraction datasets (figure[2]). When Ru was 
allowed to substitute on the A-site, the refinement converged with 7.0(3) mol% Ru 
and a slight improvement in the fit. A similar trial refinement of Nd on the B-site 
did not improve the fit. The stability of all refined models (B-site disorder, A-site 
disorder, no antisite disorder) was verified by perturbing other parameters. Because 
allowing A-site disorder improves the fit only marginally, it is important to consider 
whether the improved fit is significant, or if the improvement is merely because more 
parameters are introduced. Use of a Hamilton test [ 31i] shows the difference is 
statistically significant at <0.5% confidence interval. Indeed, with only one additional 
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Figure 4. Analysis of the synchrotron X-ray PDF reveals Nd2Ru207 is well-described 
by an average long-range structural model. 

parameter between the two models, the large number of independent measurements 
makes virtually any improvement in i?„p statistically significant. Visual inspection 
reveals only minor changes between the models, though in the limited areas where 
antisite disorder causes the largest changes, A-site disorder improves the fit to the 
experimental peak-shape (figure|3]) . The same results are obtained in samples made 
by USP and by solid state reaction. With these considerations in mind, the presence 
of A-site disorder is suggested. 

The X-ray pair distribution function (PDF) obtained by total scattering agrees well 
with the model generated by the average, long-range structure (figureQ . Refinement 
of partial Ru substitution on the A-site was attempted, but the refined occupancy 
converged to unphysical values, potentially due to strong correlation with the scale 
factor. Ideal cation site ordering Nd2Ru207 yielded a refinement R = 10.05% and 
allowing 7.0% Ru to occupy the Nd site, as suggested from Rietveld refinement, 
yielded R = 10.02%. The numerical improvement of the fit is marginal, and visual 
inspection reveals the difference between the models is much less than the level of 
noise present in the fit to the data. Rietveld refinement of Bragg scattering suggests 
partial antisite disorder, whereas PDF analysis of the total scattering shows no strong 
preference between full ordering or partial antisite disorder This likely occurs because 
the PDF refinement is strongly weighted by near-neighbor correlations, which has 
been previously noted in the system La4LiAu08 [l32ll . 

3.2. Magnetic and electrical transport behaviour 

Zero-field cooled (ZFC) and field cooled (FC) measurements of the magnetic 
susceptibility show magnetic ordering of Nd2Ru207 at Tn = 143 K (figurejs]). The 
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Figure 5. (a) Antiferromagnetic ordering of Ru''+ 4d spins leads to a small cusp 
in the ZFC magnetic susceptibility at 143 K. Nd'^+ 4f spins remain paramagnetic, 
and cause the upturn in susceptibility at lower temperatures. The separation of the 
ZFC-FC susceptibility is tentatively attributed to spin canting, which leads to weak 
ferromagnetism. (b) Scaled inverse susceptibility as a function of scaled temperature, 
as described in the text. The dashed line represents ideal Curie-Weiss behaviour, and 
the negative deviation in Nd2Ru207 is due to short-range ferromagnetic interactions. 
The Curie-Weiss fit of the high-temperature data reveals a negative Weiss temperature, 
indicating that the dominant magnetic interactions are antiferromagnetic. 




Figure 6. (a) The antiferromagnetic ordering seen in the DC magnetic susceptibility 
at 143 K is clearly visible in the xT plot, and corresponds closely with (b) the single 
anomaly in the heat capacity with a maximum at 142 K. The increase in heat capacity 
with decreasing temperature at T < 5 K is attributed to ordering of the Nd'^^ 4f spins. 
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higher-temperature region (340 K to 380 K) of the susceptibiUty data was fit to the 
Curie-Weiss (CW) equation, x = C/{T — 6cw)- The effective moment was extracted 
using the relationship /i^ff = SCk^/N, while the estimated spin-only and unquenched 
moments of Nd2Ru207 were calculated using /i^ff = 2/i|y + 2yU^j. The determined 
effective paramagnetic moment was /igff = 6.36 jj^ per Nd2Ru207 formula unit and 
the Weiss temperature was 6cw = —168 if. The /leff is close to the calculated spin- 
only value of fis = 6.8 |Ub, and significantly less than the calculated unquenched Hl+s 
of 11.9 /iB- Curie-Weiss analysis reveals a negative ©cw, indicating that the dominant 
magnetic interactions are antiferromagnetic. However, it is apparent from the first 
derivative (not shown) that the inverse susceptibility remains mildly non-linear as a 
function of temperature, indicating that Nd2Ru207 does not display pure Curie-Weiss 
paramagnetism at these higher temperatures. Consequently, the determined Weiss 
temperature and fieff should be treated only qualitatively. 

Rearranging the Curie-Weiss equation allows the scaled inverse susceptibility 
C/(x|0cw|) — 1 to be displayed as a function of T/|6cw| and provides a convenient 
way to visualize deviations from ideal Curie-Weiss paramagnetism [figure [5](b)] [l33ll . 
These deviations from Curie-Weiss behaviour are due to short-range interactions, 
and the negative deviation seen in Nd2Ru207 arises from uncompensated moments. 
Additionally, the plot provides a convenient method to visualize magnetic frustration 
(the frustration index / = 0cw/2n) [|34i] . Moderately frustrated systems tend to have 
/ > 3 [|34i1 , whereas for Nd2Ru207 / = 1.2, indicating the antiferromagnetic ordering 
is not strongly frustrated. 

Previous studies have shown that the transition at 143 K is the result of 
antiferromagnetic ordering of Ru^+ 4d spins, as an analogous transition is observed in 
Y2RU2O7, where there are no / electrons at the A site []6l [TOll . The negative Weiss 
temperature indicates the dominant magnetic interactions are antiferromagnetic; 
however, the Nd'^+ 4f spins remain paramagnetic, and cause the upturn in 
susceptibility at lower temperatures. Additionally, the history-dependence of the ZFC 
and FC measurements and the increased susceptibility at lower applied fields suggest 
there are weak uncompensated (i.e. ferromagnetic) moments. Owing to the dominant 
paramagnetic component of the susceptibility and the normalization over the applied 
magnetic field, the weak ferromagnetic interaction saturates at lower applied fields 
and its contribution appears smaller at higher applied fields in the field-normalized 
susceptibility [figure [5](a)]. Meanwhile, the ZFC-FC bifurcation of the magnetic 
susceptibility is consistent with weak ferromagnetism arising from spin-canting, and 
is described later in further detail. Below 10 K, there is another transition due to 
antiferromagnetic ordering that can be clearly seen in the xT plot [figure[6](a)]. There 
is also a corresponding increase in the specific heat that is attributed to the tail of a 
Schottky anomaly, likely resulting from the ordering of the Nd 4f spins f35\ |36l |37ll . 
The rare-earth 4f spins are expected to interact much less strongly than Ru 4d spins 
because the 4f-orbitals are contracted [l38ll . Consequently, Nd 4f spins are decoupled 
from the Ru 4d spins and order at much lower temperature [i9J . 
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Figure 7. (a) Low-temperature electrical resistivity follows a Mott 3D variable-range 
hopping model, but deviates at the magnetic ordering temperature (Tn — 143 K). 
Magnetic ordering of the Ru 4d conduction electrons leads to a change in the electrical 
transport, as shown in the derivative (b) . 



The antiferromagnetic ordering of Ru 4d spins at = 143 K causes notable 
effects in other measurements as well. The specific heat contains a corresponding 
A-type anomaly with a maximum at 142 K [figure[6](b)], consistent with a second- 
order phase transition. Additionally, the electrical resistivity near the magnetic 
ordering temperature displays anomalous behaviour. An Arrhenius-style plot shows 
a change in slope at the magnetic ordering temperature (figure[7]). Electrical 
conduction in y42Ru207 materials involves Ru 4d states, so it should come as no 
surprise that magnetic ordering of Ru 4d spins causes a marked change in electrical 
transport due to changes in scattering, as shown in the derivative of figure[7} Low- 
temperature electrical resistivity follows a 3D variable-range hopping model with p{T) 
= poexp(^)3 [|39ll40ll . with a change in the hopping barrier at the magnetic ordering 
temperature (Tn = 143 K). 

There has been considerable confusion in the literature about the nature of the 
ordering of the Ru 4d spins in Nd2Ru207 and analogous Ru-pyrochlores, in part due to 
the unusual field-dependent hysteresis present between the ZFC and FC susceptibility 
measurements, and also due to the many types of disorder and exotic phenomena that 
sometimes accompany geometric frustration in magnetic pyrochlores. In particular, 
several reports have mentioned that the glassy nature of the ordering of Ru 4d spins 
is evident in bulk magnetic susceptibility measurements Il8l |9l HOl [TTl [T2ll . However, 
neutron diffraction experiments below the ordering temperature show that Ru spins 
order with a long correlation length [|9l [TOll . This is incompatible with a glassy-state. 
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Figure 8. (a) The real component of the AC susceptibility shows a small peak 
due to the antiferromagnetic ordering of Ru 4d spins, (b) Background-subtracted 
data show the transition temperature is frequency-independent, which indicates the 
antiferromagnetic ordering of Ru 4d spins is not glassy, and that Nd2Ru207 is not a 
spin-glass. 

where there is no long-range magnetic order. This is not the only signature of a spin- 
glass that is violated by experimental evidence. Notably the temperature-dependence 
of the specific heat should vary smoothly near the magnetic ordering temperature 
[RTII . a stark contrast with the A-type anomaly at 143 K, shown in figure[6](b) . 

AC magnetic susceptibility measurements were also conducted, as another 
signature of a spin-glass is a magnetic ordering temperature that varies with the 
frequency of the applied magnetic field [l4Tl [42|| . The real, in-phase component of 
the AC susceptibility x', is shown in figure|8][a), while the imaginary, out of phase 
component (x") of the complex susceptibility was below the instrumental detection 
limit, and is thus not presented. Owing to the paramagnetic response of the Nd^+ 
that dominates the signal, the paramagnetic background was fit to the Curie-Weiss 
equation and was subtracted to yield the antiferromagnetic component of the signal 
[figure |8](b)]. The antiferromagnetic ordering temperature associated with the Ru 
4d spins is independent of frequency, which is inconsistent with canonical spin-glass 
behaviour PHI El. 

A spin-glass, by definition, must lack long-range order and display a frequency- 
dependent peak in the susceptibility PTl |42|| : neither condition is satisfied in 
Nd2Ru207. Additionally, given that the rare-earth 4f-spins are decoupled from the 
Ru 4d-spins, the nature of the magnetic ordering is not expected to differ in the 
y42Ru207 (A = Pr, ... Lu, Y) series as the identity of the rare-earth atom changes, so 
it is also unlikely that other members are spin-glasses, contrary to previous reports 
[]8l[9l[IOl[ni[l2]]. 

Adding to the confusion regarding the magnetic behaviour of y42Ru207 {A = 
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Pr, ...Lu, Y), previous studies of Nd2Ru207 have claimed there is ferromagnetic 
ordering at 20 K, and reported a hysteresis in the low-temperature field-dependent 
magnetization [7], neither of which were observed here in samples of >99% purity. 
Additionally, a later study of the specific heat of Nd2Ru207 showed multiple anomalies 
at 130 K and 20 K, in addition to the anomaly observed here at 142 K [8 ] . However, the 
ferromagnetic ordering at 20 K and the other specific heat anomalies are features of 
NdaRuOy, a secondary phase that is easily formed during preparation of Nd2Ru207 and 
to which we attribute the ferromagnetic ordering and other specific heat anomalies. 
Extensive study of NdaRuOy by neutron diffraction, magnetization, and specific heat 
measurements has shown the material undergoes a transition with a peak at 19K 
due to ordering of both Ru^+ and Nd^+ spins [|43j| . Additionally, there is a peak in 
the specific heat at 130 K, corresponding to a first-order structural phase transition. 
Although preparation of yl2Ru207 (A = Pr, . . . Lu, Y) by ceramic reaction appears 
straightforward, care must be taken at high temperatures in air to prevent slow 
decomposition of the product. To demonstrate this, Nd2Ru207 was annealed in air 
at 1040°C and 1060°C for two weeks, and led to the formation of Nd3Ru07 as a 
dominant, or single phase. In the literature, preparation of Nd2Ru207 appears to 
consistently lead to appreciable amounts of Nd3Ru07 as a secondary phase, with some 
studies showing up to 10mol% flPJ. Due to the similar properties and chemistry of 
analogous y42Ru207 systems, it is not surprising to see similar features near 20 K in the 
magnetic susceptibility and the specific heat (e.g., ferromagnetic ordering and a A-type 
anomaly) [8], as these features may originate from an yl3Ru07 secondary phase. 

3.3. Electrical transport in Ru pyrochlores 

The nature of the electrical transport behaviour in ruthenium pyrochlores (A2RU2O7) 
has been of considerable interest, as the nature of the A-site ion dictates whether 
the material will be insulating (e.g., A = Pr, . . .Lu, Y) or metallic (e.g., A = Tl, Pb, 
Bi) []3l [291 HH 051 BSl . The most intuitive model has sought to explain the different 
behaviour strictly with the A-site ionic radius, as increasing the A-site radius increases 
the Ru-O-Ru bond angle. Electrical conduction takes place in the RU2O6 network 
through overlap of Ru 4d and O 2p orbitals, so a larger Ru-O-Ru bond angle increases 
this overlap, the bandwidth, and the electrical conductivity flU |3l |45|| . 

Unfortunately, testing this relationship is not straightforward, as no lanthanide 
leads to metallic behaviour. Although the La^+ radius is similar to Bi'^+ and should 
thus be large enough to induce metallic behaviour, La2Ru207 is outside the pyrochlore 
stability-field [HI [6l [47 1 . Several studies interested in the electrical conductivity of 
these systems have examined solid solutions with Bi or Pb to overcome this hurdle 
and increase the average ionic radius to the point at which the system becomes 
metallic [[3l [22l [25l [261 [HI [48]| . However, it is not appropriate to compare rare- 
earth Ru pyrochlores with those containing Pb or Tl (e.g. Pb2Ru206.5, Tl2Ru207_y), as 
these contain either A-site ions of different formal charge in the case of Pb^+, or have 
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Figure 9. Room-temperature electrical resistivity of A2-a;Bij;Ru207 (A = Pr, . . . Lu, Y) 
solid solutions decrease smoothly with increasing average cation radius. Introduction 
of Bi drastically changes the behaviour (circles). Dotted lines are guides to the eye, 
and the thin dashed line at p = 0.01 fl cm represents the Mott minimum metallic 
conductivity at room temperature. Ionic radii of 8-coordinate 3+ cations were 
taken from Shannon 1149 II . Values of electrical resistivity were taken from references 
Ifn[3ll22ll26l[50l[5ni52ll53]|. 



contributions from overlap between empty 6s states and filled states, as in the case of 
TP+. Even an isovalent substitution of Ln^+ by Bi^+ may not allow straightforward 
comparison; structural studies of Bi2Ru207 have shown that Bi^+ atoms are off- 
centered due to the stereochemically-active 6s^ lone-pair ||5l [24ll . Indeed, Kennedy 
and Vogt, among others, point out that, while several structural or physical parameters 
of insulating A2RU2O7 (A = Pr, ... Lu, Y) pyrochlores tend to follow a simple trend 
(e.g., a linear variation in oxygen positional parameter and Ru-0 bond distance as 
a function of lattice parameter), the behaviour of metallic ruthenium pyrochlores is 
aberrant [l2l[25ll. 

To examine the influence of the A-site ionic radius on the electrical transport, 
measurements of the electrical resistivity at 300 K for several Ru pyrochlore systems 
have been gathered and are presented in figure[9| Although examining the 
relationship between the Ru-O-Ru bond angle and electrical resistivity might be more 
direct, the lack of reliable structural data makes the number of available systems less 
informative. As mentioned earlier in detail, Kennedy and Vogt demonstrated the use 
of lab XRD has led to inaccurate determination of the O atomic position Il25ll . On the 
other hand, the ionic radius is independent of the reported crystallographic data, so 
it is used as the abscissa in this case. Additionally, when examining these systems 
and structure-property relationships, it is common to show parameters plotted V5. 
lattice parameter However, the variation in lattice parameter is strongly influenced 
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by changes in the nature of bonding, as in the case of metaUic Bi2Ru207. For example, 
even though 8-coordinate Bi^+ has a larger ionic radius than Nd^+, (rBj3+ = 1.17A, 
r^^3+ = 1.109 A) [i49i], BisRuaOy has a smaller unit cell than Nd2Ru207 [l22ll . 

Casual examination of figure|9] suggests a monotonic decrease in the room- 
temperature resistivity as the ionic radius of the A site increases, when only the rare- 
earth-containing pyrochlores (A2RU2O7, A = Pr, ... Lu, Y) are considered. However, 
substitution with Bi dramatically changes the behaviour The difference in electrical 
resistivity is truly striking, and suggests that there is more than the effect of ionic 
radius when Bi is incorporated into the material. We argue that local distortion of Bi 
centres in Ru pyrochlores is responsible for the metallic behaviour and the distinct 
behaviour of Bi2Ru207 and A2-x^^x^^20-j when contrasted with Bi-free samples. 
Other researchers have pointed to the presence of additional oxygen vacancies (other 
than the ordered vacancy dictated by the pyrochlore structure) and decreased Ru-0 
bond distance that accompany metallic behaviour HSl [29ll , though these are a result, 
rather than the cause, of metallic bonding. 



3.4. High-temperature thermoelectric properties 

The electrical transport properties of rare-earth y42Ru207 (A = Pr ... Lu, Y) members 
are near the onset of metallic behaviour. This unique position is a good place to 
examine thermoelectric properties, where there is a balance between a high Seebeck 
coefficient and low electrical resistivity. Electrical resistivity and Seebeck coefficient 



of Nd2Ru207 are presented from 300 K to 900 K (figure[10|) . The thermoelectric 
properties of Nd2Ru207 do not change over three measurement cycles, despite the 
low oxygen partial pressure and high temperature during measurements. The 
room-temperature electrical resistivity is 2.0i7cm, and decreases with increasing 
temperature, as expected for a non-metal. The room-temperature Seebeck is 
promising (^220/iV/K), but quickly decreases with increasing temperature, and 
K,20^V/K above 420 K. nZT is also presented [figure[To](b)], as a proxy for the 
thermoelectric figure of merit, and is several orders of magnitude too small for 
Nd2Ru207 to be a competitive thermoelectric material. The thermoelectric figure of 
merit, ZT, is given by ZT = S'^T/{pn) and is a function of the Seebeck coefficient 
S, electrical resistivity p, thermal conductivity k, and temperature T. Many materials 
exhibit a thermal conductivity between IWm^K^ and lOWm^K^, so kZT is a 
useful proxy to compare the electrical performance of thermoelectric materials and 
estimate ZT within an order of magnitude. 

Following our investigation of Nd2Ru207 and learning of its low nZT at 
high temperatures, we realized that a simple analytical method for estimating 
the competitiveness of a thermoelectric at high temperatures would be useful for 
identifying new candidate materials. Such a method would be especially helpful 
if it did not require high-temperature measurements, as this would save time 
and resources. Many materials are currently being screened for high-temperature 
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Figure 10. (a-b) High-temperature thermoelectric properties demonstrate that 
Nd2Ru207 has poor thermoelectric performance from 300 K to 900 K. Three heating 
and cooling cycles were performed to ensure the stability of the sample; the first cycle 
is represented by hollow symbols, subsequent cycles by filled symbols, (c) Comparison 
with high-temperature p-type thermoelectrics reveals that high-performance materials 
are clustered in one region of the plot, while the performance of Nd2Ru207 is limited 
by high resistivity. This method of visualization can be used as a rapid screening 
tool to determine whether a material is likely to have reasonable thermoelectric 
performance. The radius of the circle is nZT at 700 K. Since kZT of Nd2Ru207 is very 
small, it is represented by x. Values are presented for Ca2Co205 [|54L Sio.8Geo.2 115511 . 
(AgSbTe2)o.i5(GeTe)o.85 (TAGS-85) [Ha, and Nao.o2PbTe llSTll. 



thermoelectric performance, but measurements at elevated temperatures require 
specialized instrumentation that is not widely available. Room-temperature Seebeck 
measurements can be performed quickly and with less sophisticated equipment. 
Additionally low-temperature (i.e., T<300K) Seebeck and electrical resistivity data 
are available in the literature for many materials, and could be used to quickly 
eliminate materials that are likely to have poor thermoelectric performance. 

Towards this end, a new type of plot is presented in figure[To|(c), where the 
room-temperature Seebeck (5*300 a') is plotted versus the room-temperature electrical 
resistivity (psooa)- The radius of a data point represents the magnitude of kZT at 
an arbitrary temperature that is common among the data points, and allows easy 
comparison of the expected thermoelectric performance at that temperature. When 
examining Nd2Ru207 and its relation to other high-temperature p-type materials 
at 700 K using the aforementioned survey plot [figure[TO|(c)], high-performance 
materials are clustered in one region of the plot (filled area) while Nd2Ru207 is 
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isolated. In particular, the location of Nd2Ru207 indicates the room-temperature 
electrical resistivity is too high. Because the kZT of Nd2Ru207 is very small and leads 
to a vanishingly small point, it is represented by the symbol x in figure[To|(c). 

This type of analysis is particularly useful in establishing that the performance of 
a particular class of materials may be ineffectual due to a key property being outside 
the useful range. The choice of axes is similar to a Jonker plot, in which the Seebeck 
coefficient is plotted V5. electrical conductivity. Jonker plots are used to examine the 
effect of changing the carrier concentration of a single material rather than looking 
at a field of candidate thermoelectics [l58l [59ll . Also, while Zhu et al. have shown 
that Jonker plots could be used to to estimate the peak thermoelectric power factor 
(S^p^) of a material [60], other analyses are involved, and the approach is distinct 
from the one employed here. 

4. Conclusion 

Pyrochlore Nd2Ru207 has been prepared and examined using a combination of 
structural, magnetic, and electrical and thermal transport studies. Some substitutional 
disorder on the A-site is proposed from the structural studies, but is not anticipated 
to strongly influence the physical properties. The magnetic behaviour of Nd2Ru207 
has been clarified through a combination of DC and AC magnetic measurements, and 
heat capacity studies. Despite the potential for geometric frustration of magnetism 
in the pyrochlore structure -type, we find no such exotic behaviour or glassiness in 
Nd2Ru207 and instead tentatively suggest it is a canted antiferromagnet that displays 
weak ferromagnetism. When the electrical transport properties are regarded in 
light of published data on rare-earth substituted y42Ru207 pyrochlores, it is clear 
that ionic radius plays a key role in determining electrical behaviour. However, 
the metallic electrical properties that accompany incorporation of Bi^+ on the A- 
site lies outside this description, and we suggest Bi'^+ off-centering may drive 
this anomalous behaviour High temperature measurements of the thermoelectric 
properties indicate that Nd2Ru207 has excessively high electrical resistivity for it to 
be a useful thermoelectric, despite displaying a promising Seebeck coefficient at room 
temperature. We propose a modified version of the Jonker plot as a powerful tool 
to screen candidate thermoelectric materials. We find that it is particularly useful in 
establishing that a particular materials class may be ineffectual due to a key property 
being outside the useful range. 
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